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In this work we investigate the scahng properties of quasiparticles of Pnictide with "half-Dirac" node under 
magnetic field in vortex state. By computing the density of states, we aim to find in vortex state the form of 
non- Simon-Lee scahng predicted for such system by several recent works in non- vortex state. We find by exact 
^ diagonalization of the Bogoliubov-de Genne Hamiltonian and finite size scaling a N{E) ~ ^/E power law in 
^ the case without magnetic field which agrees with analytical prediction. We consider the vortex state by first 
^ studying the hypothetical situation of uniform magnetic field without vortices and then we properly treat the 
O magnetic field-induced vortex lattice by expressing the BdG Hamiltonian in terms of superfluid velocity and 
^ Berry's gauge fields. The two calculations are shown to agree with each other. We then analyze quantitatively, 
Q the effects of anisotropic dispersion to the scaling properties of vortex lattice and show that the vortex lattice 
^ spacings will scale as dy ^ s~^d'^ where \ < rj < 1 and s ^ as compared to 77 = ^ from dimensional 
,— (Scaling analysis of non-vortex state. A very crucial prediction is also made on an upper bound to the value of 
^ 'anomalous dimension' 5 of density of states scaling with magnetic field which we find to be 5 < ^, a quantity 
Q that could not be determined conclusively by previous purely analytical works and a quantity that can be 
^ measured experimentally. 
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I. Introduction 

The discovery of the new family of high Tc super- 
conductors in iron-based compounds [1] has spurred 
tremendous theoretical and experimental activities to 
understand the mechanism of superconductivity and 
its properties. One of the main issues in this case is 
the pairing gap symmetry. It is believed that the sys- 
tem has s± pairing symmetry [2] with opposite sign 
between 5-wave gap in the electron pocket and that in 
hole pocket. However, it has been suggested from first 
principle calculations [3] that electron pocket gap func- 
tion is anisotropic and can have gapless points rather 
than being just simply isotropic 5- wave. The existence 
of gapless point has crucial importance as it can give 
rise to quasiparticles that have significant implications 
on the thermodynamic and electronic properties of the 
superconductor which can be measured in experiment. 
It is therefore of particular interest to study the impli- 
cations of this electron pocket gap anisotropy and the 
associated quasiparticles. In a ground-breaking work, 
V. Stanev et al. [4] showed that, for a critical case 
where the anisotropic gap function has gapless point 
which just touches Fermi surface, the Pnictide super- 
conductor is expected to have nontrivial scaling prop- 
erties that is distinct from the standard Simon-Lee 
scaling [5] applicable to rf-wave superconductor quasi- 
particle. This prediction was soon verified by R. Fer- 
nandes et al. [6] who extended the result to those for 
more general geometries of electron Fermi pocket and 
studied in more detail the nature of quantum phase 
transition which was also predicted in Ref. 4. The 
Simon-Lee scaling law says that for nodal quasipar- 
ticle spectrum of d-wdcve superconductor under mag- 
netic field H, the energy eigenvalue scales as En{H) ^ 
H andm the density of states scales as N{E^H) ^ 
^/H. This is true for nodal quasiparticles of (i-wave 
Cuprate superconductors, where the energy spectrum 
is isotropic around the node, that is, energy varies lin- 
early with momentum away from the node in all di- 
rections. In Pnictide, in a two-band model with ex- 
tended s± wave pairing symmetry, the electron pocket 
Fermi surface has anisotropic gap function that can be 
modeled as Ag = Aq + A'co5(40) where 9 is the an- 
gle around the M = (vr, vr) point of electron pocket in 
reduced Brillouin zone, corresponding to 2 Fe atoms 
per unit cell of FeAs lattice. With this gap anisotropy 
of electron pocket Fermi surface, we can have line of 
quantum phase transition by varying the parameter 
ratio This line of QPT occurs at ^ = 1 where 



there is transition from spectrum with node for ^ > 1 

to spectrum with zero right at the QPT line ^ = 1 to 

nodeless (fully gapped) spectrum for ^ < L One im- 
portant finding [4] is that, around the zero, occurring 
for example at momentum q = {pf^^) measured from 
the M = (tt, tt) point, the Pnictide BdG Hamiltonian 
without magnetic field 

(1) 

where A = leads to anisotropic energy disper- 

Pp 

sion, E = ^ {vfPxY + (^^§^)^. This anisotropy leads 

V Pp 
to an entirely new scaling behavior that is distinct 

from that of Simon and Lee. Following the same pro- 
cedure used by Simon and Lee, it is shown in [4] that 
in the absence of magnetic field, the density of states 
(DOS) scales as N{E) - \/E rather than N{E) ^ Ein 
Simon-Lee scaling law. The corresponding new scaling 
laws for other thermodynamics properties, for Pnic- 
tide still in the absence of magnetic field, follow ac- 
cordingly. In the presence of magnetic field applied 
perpendicular to the plane of superconductor, the sit- 
uation is much less trivial since there will be vortices. 
Authors in Ref. 4 showed that even without solving 
the vortex state Bogoliubov-de Genne(BdG) problem 
exactly, we can readily get the scaling behavior of spec- 
trum and thermodynamic properties. It is shown that 
the energy eigenvalues and density of states at low en- 
ergy (if magnetic field can create quantum states at 

2 

low energy) at low field scales as En{H 0) Hs 
and N{E = 0,i7 ^ 0) - H^^^ with a fascinat- 
ing possibility for the existence of nonzero anoma- 
lous scaling dimension S. This possibility of anoma- 
lous dimension can supposedly be associated with the 
fact that the spectrum is anisotropic. The authors in 
Ref. 4 have used at least two methods to compute the 
value of 5, one using the semiclassical argument due 
to G. E. Volovik [7], which gives ^ = — y2 another 
method based on the analogy of this problem with the 
Schrodinger problem of particle inside a box of dimen- 
sion equal to magnetic length I = which gives 
5 = |. There is apparently disagreement on the value 
of this anomalous dimension obtained from the two 
approaches. It is our goal to use another method to 
calculate this anomalous scaling dimension. So far the 



case where we have vortices in the intermediate state 
Hci < H < Hc2 has not been studied and this is what 
we pursue in this work. We here analyze the problem 
systematically by first considering a simpler though 
unphysical situation of Pnictide under magnetic field 
but without vortices in Section 11. We then do the ac- 
tual vortex state calculation in Section III and discuss 
the main results in Section IV. We especially focus on 
the computation of density of states in this paper as it 
suffices to verify the novel non-Simon-Lee scaling and 
leave the calculation of other properties, e.g. thermo- 
dynamic or magnetic properties, for future work. 



II. Density of States of Pnictide half-Dirac 
Nodal Quasiparticles under Magnetic Field 
but without Vortices 

We consider here a simpler problem of Pnictide under 
magnetic field but without vortices. We first add mag- 
netic field to the BdG Hamiltonian (1) approximated 
to the lowest order around the zero at q = (pf^O) and 
choose Landau gauge describing a uniform magnetic 
field given as Ax = Hy and Ay = to obtain the 
following Hamiltonian, 



H 



A 



A 

-vf{Px - 



c > 



(2) 



with A 



^^2^^ which represents the peculiar half- 
Dirac node character of Pnictide electron gap func- 
tion we are considering here. A different choice for 
the gauge, e.g. symmetric gauge A = — x), will 
not change the physics and the final conclusions since 
the system described by Hamiltonian (2) is gauge in- 
variant. We diagonalize this Hamiltonian exactly in 
plane wave basis on a finite size system with peri- 
odic Born- von Karman boundary condition where the 
wave vectors are given by fc^ = with L as the 
finite system size. In Fig. 1, we present the density of 
states N{E) and cumulative density of states defined 
as K{E) = /(f N{E')dE' of the half-Dirac nodal quasi- 
particles. If at low energy the density of states N{E) 
has power law dependence on E as N{E) ^ E^^ then 
the cumulative density of states has power law scal- 
ing as K{E) = N{E')dE' - E"^^^. In the case 
with magnetic field, this generalizes to K{E^ H) — 
N{E\ H)dE\ where N{E,H) is the corresponding 
density of states in magnetic field H. 
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FIG. 1 The density of states vs. energy for a Pnictide with 
half-Dirac node under zero magnetic field = (a) and 
finite magnetic field H = 2.0 but without vortices (b), ob- 
tained from exact diagonalization of Hamiltonian in equa- 
tion (2). (c) and (d) are the corresponding cumulative den- 
sity of states. 

The parameters used in producing the results shown in 
Fig. 1 are pF = 2.0, Ep = p^/2m = 2.0, Aq = 0.0001 
where is the Fermi momentum, Ep is the Fermi 
energy, and Aq is the gap amplitude. The unit sys- 
tem is chosen in such a way that one unit of mass mg 
equals mo = 10~^-'^A:^ m/s, one unit of length /q equals 
/q = 10~^m = 1/xm, one unit of velocity V{) (which fixes 
one unit of time) equals — lO^m/s which is chosen 
to be of the same order as typical electron Fermi veloc- 
ity in metals and from these follows one unit of energy 
E^ which we define to be = ^ = 0.3125ey and 
one unit of magnetic field i7o which equals i^o = ^T. 
We observe from Fig. 2 that at zero magnetic field, 
the density of states increases from zero at zero en- 
ergy while the one with finite magnetic field increases 
from some finite value at zero energy. Performing a 
finite size scaling where the result from finite size sys- 
tem is extrapolated to the infinite size limit as shown 
in Fig. 2(a), we firmly establish that in the absence of 
magnetic field, with the Hamiltonian as given by equa- 
tion (2), the density of states increases as power law 
N{E) ^ E^ where c is numerically convergent to 0.5. 
This agrees with the result of theory of non-Simon-Lee 
scaling [4,6] for the quasiparticles of Pnictide with half- 
Dirac node which predicts that for the case of super- 
conductor without magnetic field, the density of states 
should increase with energy as N{E) ^ This is 

a faster increase than that of the d-wave Cuprate case 
and this originates from the parabolic part of the dis- 
persion which turns out to give rise to more quasipar- 
ticles to appear in the vicinity of half-Dirac node than 



a full-Dirac node does in d-wave Cuprate superconduc- 
tor. 



N{E, H) = AH^{E^H-^ + B) (3) 



(a) 



where A and B are the constants with appropriate di- 
mension to be determined from numerics. This imphes 
f{E,H) = A{E^H-^ + B). The reason for this choice 
is that N{E, H) correctly reduces to N{E, H) ^ E2 at 
i7 = and to N{E, H) ^ ^ i7 at low energy 
as required from dimensional analysis. In Fig. 3, we 
show the fitting of the numerical results for N{E,H) 
with power law scaling function. 




Pnictide with half-Dirac node in the absence of magnetic Magnetic Field 

field (zero field) (b) The density of states vs. magnetic field 
for a Pnictide with half-Dirac node under magnetic field 

but without vortices, obtained from exact diagonalization FIG. 3 The power law fitting of density of states as function 
of Hamiltonian in equation (2). of magnetic field without vortices 



With finite magnetic field, where the Hamiltonian is 
given by equation (2), the density of states at = 
acquires finite value, as illustrated in Fig. 1(b). We 
will verify here whether, as predicted in [4], the low en- 
ergy weak field density of states N{E = 0, H) increases 
with the field H as power law N{E = 0, i7) - 
and compute precisely the critical exponent 7. In Fig. 
2 we show N{E = 0, H) vs. H and argue that it 
represents a power law dependence. The half-Dirac 
nodal quasiparticle density of states as function of en- 
ergy and magnetic field can be written as N{E, H) = 
Hy{E,H) = Elg{E,H) where l^l + S, f{E,H) 
and g{E, H) are dimensionless universal homogeneous 
scaling functions. The exponent of H and E to low- 
est order is determined by dimensional analysis and 
elementary analysis of half-Dirac nodal quasiparticle 
anisotropic dispersion as was first done in Ref. 4. 
To obtain the actual critical exponent for power law 
dependence on magnetic field, we fit the finite zero- 
energy value of density of states as shown in Fig. 2 
with the following scaling function, 



Fig. 3 shows that the best fitting gives a power law 
scaling exponent 7 =^ + 5 0.84 which implies an 
'anomalous dimension' of 5 = 0.5066 which is rela- 
tively close to 5 = 0.50. This is to be compared 
with computation of the anomalous dimension using 
Schrodinger particle in box argument which yields 6 = 
i or 7 = ^ and using semiclassical argument ala Volovik 
which yields 5=— ^or7=|as stated in Ref. 4. 
This value of 5 = ^ is one of the main results in this 
paper that can be of particular interest to experimen- 
talists to measure directly from experiment. It is to 
be seen later that this prediction from the hypothet- 
ical case of the strongly type II Pnictide in magnetic 
field but without vortices turns out to very remarkably 
agree with the actual vortex state analysis presented 
in the next section. 



III. Density of States of Half-Dirac Nodal 
Pnictide Quasiparticles in Vortex State 



H' = 



In the actual vortex state of superconductor under 
magnetic field, we are dealing with spatially depen- 
dent order parameter as it acquires complex phase fac- 
tor whose spatial variation describes the vortex lattice spectively as v, 
configuration. The phase (/)(r) of the order parameter 
in real space A(r) = | A(r)|e*^(^) has to satisfy the vor- 
ticity condition V x V(/)(r) = 2Ti7,YliLi ^(^ ~ ^i) and 
is in general spatially varying function in both am- 
plitude and phase parts. However, for strongly type 
II Pnictide superconductors, where inter vortex sep- 
aration is much larger than vortex core size, we can 
assume that the amplitude is constant and we are left 
with spatially varying phase function (/)(r). The Hamil- 
tonian will also contain gauge field which can be fixed 
by a particular gauge choice, e.g. Landau gauge. We 
then end up with a Schrodinger problem which is not 
exactly solvable in general, once vortices are present. 
The magnetic field is not uniform and vanishes every- 
where except at the centers (the cores) of the vortices 
where the magnetic flux penetrates through the super- 
conductor and the order parameter vanishes. We treat 
this periodic arrangement of vortices by rewriting the 
BdG Hamiltonian in terms of superfluid velocity and 
the Berry's phase gauge fields and perform a singu- 
lar gauge transformation to deal with multi-valuedness 
of the wavefunction due to tt phase difference (minus 
sign) it acquires as the quasiparticle winds around a 
vortex line[8,9,10]. From Hamiltonian in equation (2) 
but now in an actual vortex state, 



Vf{Px+CLx) -\-mvFVsx 
{Py{Py-\-2mVsy-\-Say)) 



8Ao 



-vf{Px + ttx) + rnvpVsx 
(7) 



where the superflow and gauge fields are defined re- 



s - ^(v^ + vf) anda- 

|(v</.^ - s/ct>B) with = ^{nvct>^ 

write H' in terms of Pauli matrices, 



m 



2 \ s 
c / 



We can 



8Ao^ 



= mvFVsx(^o+yF{Px+cix)(^3-\ — ^ {Py{Py+'2rnVsy+3ay))ai 

Pf 

(8) 

We progress by writing the Fourier representation for 
superflow velocity field 



'(r) = 



27r ^ 



zQ.R[^+iQ.r 



(9) 



H 



2m 



A 



A 

2m 



+ ep 



(4) 



where fi = A^B. We define our vortex lattice as a 
square vortex lattice with sublattices A and B in such 
a way that = — Rq and R^ = Rq. The depen- 
dence on magnetic field H enters through lattice spac- 
ing of the vortex lattice which in general is of the order 
magnetic length d Ib = ^J^- The definition (9) 
alone already represents the periodicity of vortex lat- 
tice. The spectrum of the velocity field is thus given 
by V^(Q) - ^e-^Q-^o, the magnitude of which 
decreases with wave vector Q. From this Fourier rep- 
resentation, we can easily derive the Fourier represen- 
tation of velocity and Berry gauge fields and hence 
the Hamiltonian in Fourier space. We can separate 
the Hamiltonian (7) into two parts. The first part, the 
unperturbed one, is 



i(j){r) 



}} + ^A(r)a2</>(r) is 



where A = ^{dy,{dy,Aoe 

the gauge- invariant gap operator, we perform Franz- 
Tesanovic transformation [8,9,10] on this Hamiltonian, 
defined as 



^^0 



v 



VFPx 
8Aop| 



-VFPx 



(10) 



While the second part, the perturbed Hamiltonian, is 



H U-^HU 



(5) i^l 



VFa^ + mvFVsx 
^Py{2mvsy -\-3ay) 



VFax + mvFVsx 



where 



(11) 



(Qi4>A{r) \ We include the terms with py in the off-diagonal ele- 

Q-i^B(r) j V ) ments here because they are first order in the weak ve- 

locity and gauge fields and are thus weak with respect 
and expansion around q = (pf^O)^ we have to the 'unperturbed' Hamiltonian. Upon performing 



Fourier transformation using (9) and with the defini- IV. Scaling Analysis of Low Energy Density of 
tion of full lattice superflow velocity and gauge fields, States Dependence on Magnetic Field in 
the BdG equation in vortex state is Vortex State 



Q) 



E 



(12) 



where = Q.Rq. For Bravais lattice, as is the case 
for our square vortex lattice, Rq = + 
with Q = (^^i + ^^2) and aihj = 2Ti5ij. This implies 
that the only reciprocal lattice wave vectors that con- 
tribute to the quasiparticles dynamics in vortex state 
are given by (Qx^Qy) = 27r(^,^). Notice care- 
fully that we have to assign different variables for the 
lattice spacing in x and y directions dx and dy re- 
spectively to anticipate the anisotropy of the system 
that will manifest also in vortex state. This is a cru- 
cial difference from the d-wdive case of Cuprate where 
the X and y directions are scaled isotropically and so 
we can use the same lattice spacing d if we assume 
a square vortex lattice. This again just reflects the 
anisotropy of the system that we have.We therefore 
have 6 = 27r{^, ^).l{dx,dy) = f(m^ + m^). The 
original Schrodinger equation for this problem is then 
eventually equivalent to Bloch equation [11] which has 
the following generic form, 



/?2 

(;^(k - K)^ - ^)ck_K + E Uk'-kck-k' = (13) 



Here k is plane wave vector inside the first Brillouin 
zone satisfying Born-von Karman (periodic) bound- 
ary condition while K is the reciprocal vortex lattice 
vectors which are infinite in number if we consider 
infinite system. Bloch equation describes a system 
with periodic potential and here the periodic potential 
is given by the superflow velocity and Berry's phase 
gauge fields arising from the periodic vortex lattice. 
It needs to be pointed that such periodic potential 
representation is valid in this problem only if the aver- 
age magnetic field is zero as is the case here. Solving 
Bloch equation simply means diagonalizing the Hamil- 
tonian in momentum space expanded with respect to 
the Bloch wave vectors in the first Brillouin zone of 
vortex lattice which is accomplished here numerically. 



In this section, we will describe the effects of anisotropy 
to the quasiparticles scaling in vortex state. Looking 
at the approximated Hamiltonian in equation (2), 



H 



A 



A 

-vf{p. - ^: 



(14) 



with A = ^^^y- we see that the Hamiltonian will 
be homogeneously-scaled if we rescale the variables as 
follows, 

Px = spx 
Py = s^^^Py 
y = sy' 



X = rx 



(15) 



where r is yet undetermined and is to be set later. This 
way, the Hamiltonian can be written as i7 = s^H^ for 
some exponent n. This suggests that in all later calcu- 
lations, a proper treatment should keep the Hamilto- 
nian writable in such form and so transformations (15) 
should be used. The eigenvalues will correspondingly 
transform as = s^E^We can see here that the mo- 
mentum space is anisotropically-scaled and so is the 
real space, in general. This just reflects the anisotropy 
of the system. This has crucial importance to the vor- 
tex state calculation since it now implies that the re- 
ciprocal vortex lattice should also obey (15). That is, 
the momenta in px and py directions in reciprocal lat- 
tice are not scaled the same way and so this implies 
the vortex lattice spacings in x and y directions are 
also rescaled differently. Let us consider again the vor- 
tex state Hamiltonian (7) which basically describes the 
quasiparticles dynamics under periodic superfluid ve- 
locity and Berry's phase gauge fields to the lowest (1st) 
order in magnetic field strength. Looking at the defini- 
tion of the fields given just after equation (7), consider- 
ing Vsx^ we would naively expect that the x coordinate 
transform as x = rx' = ^x. But the off-diagonal terms 
of the vortex state Hamiltonian requires us to have for 



ay and Vgy^ ay 



ciy^ ^sy = s^^'^'^sy Order to have 



H = s^H\ Let us now consider the explicit defini- 
tion of superfluid velocity component Vgy as given by 
its Fourier series representation in equation (9) which 
can be rewritten as follows. 



27r ^ — isC 



-iQ.Rff+iQ.r 



iQ.Rj^+iQ.r 



(16) 



Focusing on the Fourier coefficient given by term in 
the fraction and using the value of reciprocal wave 
vector for a periodic boundary condition given after 
equation (2), {QxiQy) — 27r(^, the denomina- 
tor immediately suggests that the vortex lattice spac- 
ings satisfy d!y ^ s~^d^^. We aim to have a rescaling 
of the form Vsy ^ s^^^Vsy and this will require spac- 
ing d to satisfy d' ^ s^^^d while imposing condition 
(ttx, Vsx) ^ s{a'^, Vsx) on {a^, Vgx) appearing on the di- 
agonal elements will lead to another condition of the 
form dy s~^d'rj.. We see that we have two different 

scaling conditions for the same quantities; d^ ^ s~^d^^ 
and dy ^ 



s~^d^^. Apparently, to this lowest order of 
approximation for vortex state Hamiltonian, there is 
no a set of variable scaling transformations that can 
make the Hamiltonian homogeneously scaled. We in- 
terpret this as a new requirement that the actual scal- 
ing in the true vortex state should eventually take a 
form dy ^ s~^d'^ where ^ < r] < 1. This is the pre- 
cise rescaling that we have to use in doing vortex state 
calculation. This result is to be compared with that 
of dimensional scaling analysis for non-vortex state 
of Ref. 4 where redefinition of coordinate variables 
= x/d and = yd/l^ (using Ref. 4's notations, 
d = {l^/a)^^a = and / = . /-%^)implies coordi- 



nate rescaling of the form X — S X and y' ^ s'2y and 
this corresponds to 77 = |. The difference in the expo- 
nents from those of our vortex state is to be noted here. 
Thus the dimensional scaling analysis exponent from 
Ref. 4 actually corresponds to the lower bound of the 
true vortex state exponent. This point is the key dis- 
tinction between dimensional scaling analysis and true 
vortex state scaling. It is to be noted that the vortex 
lattice in real space, assuming a square lattice, is not 
physically deformed, e.g. into a rectangular lattice due 
to this anisotropic scaling. This peculiar result never- 
theless takes effect (and thus must be correctly taken 
into account) instead internally in calculations involv- 
ing the Hamiltonian such as density of states or band 
structure calculations. One interesting open problem 
is then what will happen if the vortex lattice is in- 
finitely anisotropically rescaled, i.e. the dimension in 



one direction is infinitely larger than that in the per- 
pendicular direction with regard to power law of those 
thermodynamic or electronic properties. 

We present in Fig. 4 the result of calculations for den- 
sity of states of half-Dirac nodal Pnictide quasiparti- 
cle in vortex state. In the plots that follow, we use 
vp — 100.0 and = 100.0 where vp is the Fermi ve- 
locity, = ^ is the " gap velocity" . The unit system 
is the same as that used in Section II. In the following 
plots, d represents the vortex lattice spacing which is 
taken to be of order of magnetic length. 
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FIG. 4 (a) The density of states and (b) the (low energy) 
cumulative density of states of half-Dirac nodal Pnictide 
quasiparticles in vortex state with vf = 100.0, = 100.0, 
and d = 10 (c) The critical exponent /3 defined in equation 
(18) as function of lattice spacing d 

It is found that in most range of the strengths of mag- 
netic field, the density of states in vortex case vanishes 
at zero energy despite relatively strongly fluctuation at 
finite energy. This behavior agrees perfectly with all 
past works on nodal d-wave Cuprate superconductor 
[8, 9, 10]. With such fluctuating density of states, it is 



crucial to define a correct scaling law for the density 
of states that we propose to be valid even in this vor- 
tex phase. In this case, we assume the global trend 
or average of the fluctuating density of states at any 
given energy and magnetic field to obey a power law. 
We then write the power law as follows, 

N{E,H) = aE^mE-^)^ (17) 

which perfectly agrees with the observed behavior of 
N(E^ H) shown in Fig. 4(a). The justification for this 
scaling form is that N{E^H) correctly reduces to the 
dimensional analysis result for 7 = 1. As before, we 
better consider the cumulative density of states which 
will give entirely equivalent conclusion as the density 
of states does. The cumulative density of states from 
(17) is then given by 

K{E,H) = / N{E',H)dE' = ^^E^^^'^^Hi 
Jo 3-/5 

(18) 

and in this form the critical exponent (3 is now what 
we are looking for. In Fig. 4(c) we show the expo- 
nent ;5 as a function of lattice spacing d. Since as 
mentioned before the lattice spacing of the order of 
magnetic length d ^ Ib — yj^-) we see that the crit- 
ical exponent is very small at very strong magnetic 
field. But then, as we decrease the magnetic field, the 
critical exponent increases and finally, at the limit of 
very large d Ib = y^^? or very weak magnetic field, 
the density of states critical exponent settles down to a 
stable value f3 ^ 2.5. The equilibrium value of critical 
exponent /3 ~ 2.5, looking at the form of scaling law 
in equation (17), implies a power law dependence on 
magnetic field with critical exponent 7=3^ 0.8333. 
This comes in remarkable agreement with the result 
from Section II where we considered the hypothetical 
situation of Pnictide under magnetic field but with- 
out vortices and obtained 7 ^ 0.84. One can say that 
the calculation in Section II and that in Sections III 
and IV represent two different levels of solution to the 
same vortex state problem which give correction to 
the dimensional scaling analysis and which turn out 
to give scaling law which clearly agree with each other 
between the first two. This eventually supports the 
'hand waving' argument that the authors in Ref. 4 
used to derive the proposed non-Simon-Lee scaling 
law. This striking agreement between the two val- 
ues obtained from two completely different methods 



provides yet another support for the validity of our 
result for this experimentally measureable and theo- 
retical crucial 'anomalous dimension' that accompa- 
nies the novel non-Simon-Lee scaling law for half-Dirac 
nodal quasiparticles. This result for anomalous dimen- 
sion which gives correction to the critical exponent ob- 
tained from simple dimensional analysis seems to be 
rather large, but since in this work we treat the vortex 
state properly, we can argue that the value of 'anoma- 
lous dimension' predicted here should be more valid 
than the alternative results estimated in Ref. 4. Since 
the critical exponent computed in this work converges 
to a value that is smaller in the infinite size limit than 
the value for a finite size system, as illustrated in Fig. 
2(a), it is more accurate to say that the anomalous 
dimension 5 obtained here is an upper bound to the 
actual value one might obtain in experiment while the 
lower bound is S = — ^. This result has direct rele- 
vance to experiment where this quantity can be mea- 
sured directly. However, since half-Dirac node is a crit- 
ical point that can only be obtained by fine tuning, it 
is highly dependent on the experimental feasibility to 
adjust control parameters in order to be able to verify 
the predictions made in this and previous works. 

V. Conclusions 

In conclusion, we have verified numerically the non- 
Simon-Lee scaling law predicted for Pnictide with half- 
Dirac node on electron pocket in the case without mag- 
netic field. The corresponding scaling law in the pres- 
ence of magnetic field but without vortices in terms of 
density of states has been calculated by exact diago- 
nalization of BdG Hamiltonian. The obtained power 
law directly gives the first level of correction to crit- 
ical exponent from of dimensional scaling analysis in 
terms of 'anomalous dimension'. The computed den- 
sity of states in vortex state with proper treatment of 
the superfluid velocity and Berry's phase gauge fields 
of the vortices follows a power law that also leads to 
a prediction of the value of anomalous dimension and 
serves as the next level of correction to the dimen- 
sional scaling analysis result. The difference between 
scaling anisotropy in true vortex state and that from 
dimensional scaling analysis is elucidated. In sum- 
mary, this work provides another verification and fur- 
ther progresses the investigation of the hypothesized 
non-Simon-Lee scaling proposed for such critical half- 
Dirac nodal Pnictide. 
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